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lon-beam-plasma-interaction plays an important role in the field of Warm Dense Matter (WDM) 
and Inertial Confinement Fusion (ICF). A spherical theta pinch is proposed to act as a plasma target 
in various applications including a plasma stripper cell. One key parameter for such applications is 
the free electron density. A linear dependency of this density to the amount of energy transferred 
into the plasma from an energy storage was found by C. Teske. Since the amount of stored energy 
is known, the energy transfer efficiency is a reliable parameter for the design of a spherical theta 
pinch device. The traditional two models of energy transfer efficiency are based on assumptions 
which comprise the risk of systematical errors. To obtain precise results, this paper proposes a new 
model without the necessity of any assumption to calculate the energy transfer efficiency for an 
inductively coupled plasma device. Further, a comparison of these three different models is given 
at a fixed operation voltage for the full range of working gas pressures. Due to the inappropriate 
assumptions included in the traditional models, one owns a tendency to overestimate the energy 
transfer efficiency whereas the other leads to an underestimation. Applying our new model to a 
wide spread set of operation voltages and gas pressures, an overall picture of the energy transfer 
efficiency results. 

PACS numbers: 52.55.Ez, 52.80.Tn, 52.25.Fi, 52.50.Dg 


I. INTRODUCTION 

In former works Z-pinches have been used as plasma 
targets in scientific fields like the research on Inertial 
Confinement Fusion, High Energy Density Physics and 
accelerator development Plenty of achievements 

have been obtained in these scientific fields. As the Z- 
pinch is a capacitive discharge between two electrodes, 
erosion plays an important role in the lifetime limita¬ 
tion of such devices. Therefore the Z-pinch is not the 
optimal choice for high repetition rate applications e.g. 
as a plasma stripper. However, an inductively coupled 
plasma device has the advantage to use external coils 
without contact with the plasma to initiate and maintain 
the gas ionization. The so called theta-pinch combines 
the high lifetime with the free electron density compara¬ 
ble to those of Z-pinches M- 

Gonventional small scale theta pinch devices have been 
broadly investigated in various applications [13,[HI- The 
maximum energy transfer efficiency was determined to 
be 59% [13. However, in 2008, a spherical theta pinch 
with the significantly advanced energy transfer efficiency 
up to 80% was realized by C. Teske and J. Jacoby Q. 
Later on, the efficiency was raised to 86% Q. All these 
results for the energy transfer efficiency were calculated 
out by applying the model first proposed by S. Aisenberg 
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in 1964 [13 . 

The key parameter for the plasma stripper is the free 
electron density. To achieve the higher free electron den¬ 
sity, the method is to upgrade the device. However, the 
remaining difficulty is how to forecast the electron den¬ 
sity after upgrading. Fortunately, a linear scaling law 
between the free electron density and the deposited en¬ 
ergy in plasma was shown in Q. To derive the electron 
density, the knowledge of the precise energy transfer ef¬ 
ficiency values is essential. Reviewing the former energy 
transfer efficiency models [3, [Tol - [l3 , all of them contains 
an inaccurate assumption that the reflected plasma resis¬ 
tant is constant. Such treatment is to solve the differen¬ 
tial equation which will be well described in this article. 
To obtain the precise results, a new model based on the 
experimental measurement is proposed. 


II. EXPERIMENTAL SETUP 

The basic working principle and experimental setim 
have been described already in preceding papers [7|- 
iiillll- In this article, to apply as a plasma stripper, 
the spherical theta pinch is reconstructed and upgraded 
as shown in fig. [TJ The device axis is positioned hori¬ 
zontally to match the ion beam direction. To achieve a 
higher electron density, six capacitors with each a capac¬ 
itance of 25 fjF are used, being connected parallelly in 
three groups of two capacitors connected in series. This 
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FIG. 1. Spherical theta-pinch setup for the plasma stripper 
application 
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FIG. 2. LRC equivalent circuit for the plasma stripper set-up 


setting results in a total capacitance of 37.5 fi F with a 
maximum charging voltage of 18 kV. To satisfy the de¬ 
mand of a long life-time even at high repetition rates, the 
maximum operation voltage is limited to 14 kV. Never¬ 
theless, this means a maximum of 3.7 kJ stored energy 
in the capacitors which is almost twice as much as in the 
previous setup Q . Correspondingly, the thyristor-stack 
is also improved for adapting the operation voltage of 14 
kV. Besides, the bigger glass vessel of the 6000 dis¬ 
charge volume is closely encircled by a seven-turns coil. 
The Capacitors, the switch, the transmission line and the 
coil form a typical LRC (inductance, resistance and ca¬ 
pacitance) circuit with a resonance frequency of about 10 
kHz. 

A mixture of mainly Argon and 2.2% Hydrogen is used 
as a working gas and the pressure is measured by a full 
range pressure gauge. As a diagnostic tool, a fast photo 
diode is installed to monitor the plasma. In addition, a 
Rogowski coil measures the circuit current in the trans¬ 
mission line. These two diagnostic signals are recorded 
by an oscilloscope setting at two sample rates of 1 and 
2.5 MHz. 


III. MODELS OF ENERGY TRANSFER 
EFFICIENCY 

To calculate the energy transfer efficiency of theta 
pinch devices, S. Aisenberg et al. first proposed an equiv¬ 
alent circuit as Fig. [5] [iG . [l^ . This is a typical LRC 



FIG. 3. Typical current signal shape and corresponding fit 



FIG. 4. Light signal from the plasma with assumed plasma 
start time 


circuit. In this circuit, Li and Rq are undesired parasitic 
inductances and resistances respectively. Lq is the induc¬ 
tance of the coil while represents the reflected plasma 
resistance which is the transformed value from the real 
plasma resistance into the LRC circuit. Cq is the total 
capacitance of the capacitor bank with the initial charg¬ 
ing voltage Uq. Under the frame of this equivalent circuit, 
the traditional and the new models are derived as follows. 


A. Traditional models 


1. S. Aisenberg’s model 


As Fig. [2] fulfils Kirchhoff’s voltage law, the circuit 
equation is expressed as 

dlit) , Cl(t')dt' 

( Lo + Li)-^ + (Ro + R,)m + y = 0.(1) 


To solve this differential equation, Rp is assumed to be 
constant during the whole discharge. Then Eq. © is 
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Analogously, for a free oscillation, the energy dissi¬ 
pated in the parasitic resistance equals the whole stored 
energy B: 


E = 





-IqRoTq, 


( 10 ) 


where 2 to is the decay constant for the free oscillation. 

Finally, the energy transfer efficiency rj from the ca¬ 
pacitors to the plasma is 


-q = 


= 1 _ Zp 


E 


To 


( 11 ) 


FIG. 5. Experimental value for the parasitic resistance 


2. C. Teske’s modification 


transformed to 

( Lq + Li)i[t) + [Ro + Rp)i{t) -I- —I{t) = 0, (2) 


In this modification, it is assumed that the plasma 
start time tp indicated by the light emission from the 
plasma is not the same as the current. Hence, the energy 
transferred to the plasma is 


Generally, there are three different solutions depending 
on the value of A = {Rq + i?^)^ — 4 ■ 


_Ro+^ ( yZA \ 

7(0 = 7„e ws sin , A < 0, (3) 

lit) = hte , A = 0, (4) 


lit) = I^e 2{ro+ri)* 


\/A j. — ^ 

g 2(I<o + 7-i) — g 2{I/o + -l^j) 


A > 0.(5) 


Considering the current curve shape in Fig. [31 which is 
a damped oscillation, Eq. ([3]) is the only proper solution 
which is simplified to 

lit) = sin(a;t) (6) 


f+ CO 


Wp = 


I^it)R'dt 


cos(2a;tp) 


1 + . n, . sm(2wtp) 


4a;V2 -h 1 


- 1 - 1 


( 12 ) 


Usually, the plasma ignites after the first half wave which 
suggests tp ~ tt/uj as Fig. [4] shows. Besides, ojTp is as¬ 
sumed to be much greater than I. As a result, the energy 
deposited in the plasma is approximated as 


Wp 



(13) 


with Tp = , uj = , , lo = 3 , where 

2Tp is the decay constant while w is the oscillation fre¬ 
quency of the circuit. In the case of iRo+R'p)^ <C 4 , 

ui is approximated as 


UJ = 


C'o(Ao 


( 7 ) 


Correspondingly, the modified energy transfer efficiency 
formula is attained as 



(14) 


B. New model 


Further, the energy consumed in the plasma is 

7--1-00 

Wp= / I^it)R'dt 


— 2^0^pTp I 1 


I 


407^2 + lj ' 


( 8 ) 


In the case of ujTp > 5, the approximation for Wp within 
an error of 1% is 


Wp 


^iQ-npTp. 


(9) 


In the preceding models, the expressions for the en¬ 
ergy transfer efficiency are both based on the assumption 
that the reflected plasma resistance is constant. Since 
the light signal oscillates with time shown in Fig. [U nei¬ 
ther the plasma parameters nor the value of the reflected 
plasma resistance can be constant. Hence, the traditional 
models contain the risk of deviating from the true value. 

To avoid this risk, a new model without any assump¬ 
tion is established. The key point for this new model is 
that the parasitic resistance is constant as Fig. |S] shows, 
for the experimentally measured values. Moreover, since 
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FIG. 6. Decay constant for the case of free oscillation at 6 kV 


FIG. 7. Efficiencies obtained with the three methods at 6 kV 


the capacitance and inductance do not consume the en¬ 
ergy, the total energy is completely dissipated in the par¬ 
asitic resistance and the reflected plasma resistance. As 
a result, the energy deposited in the plasma is 

OO 

Wp = E-RoAtY,lHU), (15) 

i=l 


where E is the total energy stored in the capacitors; At 
is a constant time interval between ti and ti+i; I{ti) is 
the measured current value at ti. Even though At is 
considered to be infinitesimally small, the sum is used to 
take into account the finite time resolution of the oscillo¬ 
scope. For our experiment, the substitution error is less 
than 0.01%. 

For the free oscillation, the total energy which is com¬ 
pletely dissipated in the parasitic resistance is 

OO 

E = RoAtoY,IUt^)■ ( 16 ) 

i=l 

Here, loiti) is the measured current value at ti. Ato is 
again a constant time interval between ti and U+i. 

Thus, the energy transfer efficiency is written as 


loiti)' 


(17) 


IV. RESULTS AND DISCUSSION 

In this article, an operation voltage is applied from 6 
to 14 kV. For each operation voltage, the gas pressure 
varies in the range of 0.6-160 Pa and an additional ’’zero 
shot” in the order of 10“^ Pa. This ’’zero shot” repre¬ 
sents the free oscillation of the current. The black line 
in Fig. [5] shows this kind of free oscillation at 6 kV. The 
red line displaying a good overlap with the black one is 
the fitting curve. Gained from the fitting parameters, the 
decay constant of this free oscillation is 1.02 ms. Analo¬ 
gously, the decay constants for the different pressures are 
obtained. 



FIG. 8. Energy transfer efficiencies of a complete pressure 
and voltage set 


Directly applying Eq. (HU from the S. Aisenberg’s 
model on these decay constants, the corresponding en¬ 
ergy transfer efficiencies for the different gas pressures are 
calculated. For the C. Teske’s modified model, another 
parameter besides the decay constants is the plasma start 
time tp. This is obtained from the difference between the 
beginning of the light emission in time and the current 
start time. Then, adopting Eq. (HI , the modified energy 
transfer efficiencies are also computed. Both of the pre¬ 
ceding models need the plotting and fitting procedures 
which cause propagating errors. However, for our new 
model, there is no need to plot and fit the data because 
Eq. (fTfl) is straightforwardly applied on the measured 
current values from the data sheet. Consequently, the 
energy transfer efficiencies are derived from this. 

A comparison is made among the results from these 
models at 6 kV shown in Fig. [71 Their tendencies are 
extremely similar. At low pressures, the curves rise very 
fast. Then they slow down at medium pressures. Fi¬ 
nally, all of them suddenly drop down to the point of 
60 Pa which is approximated the breakdown threshold. 
The biggest difference among the models can be found 
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at medium pressures over a wide range. The S. Aisen- 
berg model shows the highest efficiencies whereas the C. 
Teske’s gives the lowest in this pressure range. The max¬ 
imum energy transfer efficiency for the S. Aisenberg’s 
model is 78% while 52% for C. Teske’s modified model. 
In our new model, the value for the maximum energy 
transfer efficiency is 71%. Due to no assumption made 
in the new model, it is considered reliable and precise. 
Hence, the S. Aisenberg’s model overestimates the en¬ 
ergy transfer efficiency. On the contrary, the C. Teske’s 
modification model underestimates the values. 

These deviations for the traditional models are mainly 
caused by the assumption that the reflected plasma re¬ 
sistance is constant, which is untrue. The truth is the 
reflected plasma resistance varies with time as the light 
signal in Fig. |4] shows. However, this assumed constant 
reflected plasma resistance leads to an approximately av¬ 
erage value for the true time-dependent resistance. Re¬ 
garding the energy consumed in the plasma, it does not 
only depend on the resistance but also the square of the 
current. Namely, the square of the current is the weight¬ 
ing factor of the resistance for calculating the transferred 
energy. Due to the oscillation decay behavior of the cur¬ 
rent, the weighting factor for the first half wave is bigger 
than for the other half waves. Unfortunately, the resis¬ 
tance for the first half wave is almost zero indicated by 
the light signal shown in Fig. |4l Consequently, the calcu¬ 
lated energy transfer by applying the average resistance 
is overestimated. This accounts for the S. Aisenberg’s 
model’s overestimation. For C.Teske’s modified model, 
this first half wave of the current is omitted for calcu¬ 
lating the transferred energy according to the Eq. (HU. 
However, the calculation of the average resistance still 
counts in this first half wave where the reflected plasma 
resistance is almost 0. This results in the underesti¬ 
mation of the average resistance. Correspondingly, the 
transferred energy calculation is underestimated. Be¬ 
sides, the neglection of the sine term in Eq. (USD also 
contributes to this energy underestimation. For instance 
at 30 Pa and 6 kV, the maximum deviation for this sine 
term is about 3%. 

Due to the advantage of our new model and the faults 
of the traditional models, our new model is adapted 
to calculate the energy transfer efficiency from the 
measured current values. Fig. |S] shows all the measured 


energy transfer efficiencies vary with the pressures of the 
different operation voltages. In our measured pressure 
range, the breakdown threshold is only observed at 
the lower voltage of 6 and 7 kV. For each voltage, 
the pressure corresponding to the peak value of the 
efficiency is defined as the optimal pressure which shows 
the maximum transfer efficiency. Comparing all of the 
plots in Fig. |5] , it is found that the optimal gas pressure 
shifts to higher values when the voltage is increased. 
The maximum energy transfer efficiency found here is 
73% which occurs at around 80 Pa under the maximum 
voltage of 14 kV. 

V. CONCLUSION 

A new model is proposed to calculate the energy trans¬ 
fer efficiency of theta pinch in general and especially for 
the spherical theta pinch. As this new model contains no 
assumptions, the results obtained with it are considered 
exact. 

The deviations of the traditional models of S.Aisenberg 
and C.Teske from the real values are mainly caused by 
the improper assumption that the reflected plasma re¬ 
sistance is constant. The first leads to overestimation 
whereas the latter to underestimation. The correspond¬ 
ing explanations are given respectively. 

The latest setup of the spherical theta pinch that is 
used for the presented investigations has a maximum en¬ 
ergy transfer efficiency of 73% for an Argon-Hydrogen gas 
mixture which is promising especially for high repetition 
rate applications e.g. flash lamps or a plasma stripping 
device. 
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